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Abstract 

The work is devoted to the critical analysis of theoretical prediction 
and astronomical observation of GR effects, first of all, the Mercury's 
perihelion advance. In the first part, the methodological issues of ob- 
servations are discussed including a practice of observations, a method 
of recognizing the relativistic properties of the effect and recovering 
it from bulk of raw data, a parametric observational model, and fi- 
nally, methods of assessment of the effect value and statistical level of 
confidence. In the second part, the Mercury's perihelion advance and 
other theoretical problems are discussed in relationship with the GR 
physical foundations. Controversies in literature devoted to the GR 
tests are analyzed. The unified GR approach to particles and photons 
is discussed with the emphasis on the GR classical tests. Finally, the 
alternative theory of relativistic effect treatment is presented. 

Key words: Kepler's motion, General Relativity, Mercury, ad- 
vanced perihelion. 
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1 Introduction 

1.1 History in brief 

The history of the problem of Mercury perihelion advance is presented in 
PQ, [2] and elsewhere. Basically, it started with the work by the French 
astronomer Urbain Le Verrier (1859) who calculated a perturbation effect 
caused by interaction of Mercury with other planets and found the anomaly, 
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- the perihelion advance unexplained by the Newton's theory. He thought 
that the effect could be explained by planetary matter, though it was not 
confirmed later on. The American astronomer Simon Newcomb [3] (1882) 
verified the effect and assessed its value being about 39" (arc-seconds) per 
century. He also produced detailed tables of planet orbits (ephemerides) and 
suggested that a solution of the perihelion problem could be some alteration 
to the Newtonian law. Many great scientists before and after Einstein's work 
(1915) came to the scene with attempts to develop different approaches to 
the problem and improve observational techniques and calculational (per- 
turbative) algorithms, among them the American astronomer and dedicated 
researcher in the problem of Mercury perihelion advance, Gerald Clemence. 

In a long run (from Le Verrier to Einstein), there was a wide range of 
admissible values of "the anomaly" (say, from 5 to 50 arc-seconds per cen- 
tury). Probably, a strive for "a true explanation" of the anomaly prevailed 
over work on determination of "the true anomaly" value from astronomical 
observations, still of poor precision. With Einstein's GR advent, however, 
hot discussions of the problem were not ceased. Many Einstein's promi- 
nent counterparts remained highly critical about the Einstein's theory or did 
not accept it at all. When Einstein received a Nobel Prize (1922) for great 
achievements in Physics, there was a special note saying that the prize was 
given "without taking into account ... relativity and gravitational theories" 
i- 

One should bear in mind that in the historic time till, say, mid 50s, a con- 
duction of computations of planetary orbits with mechanical or electrome- 
chanical calculators was an extremely hard job. The electronic computer era 
began in about 1945, while the state-of-art maturity came after Clemence's 
activity discontinuation. In spite of technical difficulties, Clemence made 
a tremendous contribution to the problem study (early 40s till late 60s), 
especially in matching physical theories and astronomical observations and 
achieving better clarity of the problem and higher precision of data analy- 
sis. He rechecked Le Verrier's and Newcomb's results and updated planetary 
data [4J (1943), [5] (1947), [6J (1948), [Tj (1965) . It was clearly understood 
by Clemence (and pretty much so by his predecessors) that the accuracy 
of the assessment of "the anomaly" is restricted by its smallness (less than 
1/100 of "perturbation background"). There were unavoidable methodolog- 
ical problems such as a connection of time scales in different observations as 
well as determination of inertial reference frames the effect must be referred 
to. Basically, he stated that, given realistic uncertainties of input data and 
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model parameters, the planetary ephemirides can be, in principle, adjusted 
by minimization of the anomaly gap in concordance with the General Rel- 
ativity prediction. This was a solid result but still not a resolution of the 
problem because of a number of controversial issues remained unresolved. In 
1947, Clemence noted [5J (1947) : 

"According to general theory of relativity, the elliptical orbit of a planet 
referred to a Newtonian frame of reference rotates in its own plane in the 
same direction as the planet moves... The observations cannot be made 
in a Newtonian frame of reference. They are affected by the precession of 
the equinoxes, and the determination of the precessional motion is one of 
the most difficult problems of observational astronomy. It is not surprising 
that a difference of opinions could exist regarding the closeness of agreement 
of observed and theoretical motions... I am not aware that relativity is at 
present regarded by physicists as a theory that may be believed or not, at 
will. Nevertheless, it may be of some interest to present the most recent 
evidence on the degree of agreement between the observed and theoretical 
motions of the planets" . 

The situation in observational astronomy radically changed with the tech- 
nological boost started in 1950s and continued since then. After the WW2, 
astronomers armed with advanced observational tools and computer tech- 
nologies took a chance to pursue ambitious unprecedented projects in obser- 
vational astronomy and celestial mechanics, such as the international celestial 
reference frame and coordinate system, the Earth rotation service with coor- 
dinate system transformations, the database for a wide range of astronomical 
and physical constants, parameters of theoretical N-body gravitational mod- 
els, etc. 

The initiative was launched to include the GR theory in the form of the 
Parameterized Post Newtonian (PPN) approximation. It was motivated by 
needs Astronomy and Astrophysics, space missions and other national tasks 
related to time and position precision in different reference frames. Such tasks 
could be fulfilled only with the use of contemporary electronic technologies. 
In particular, it became possible to study errors in approximate solutions of 
N-body problem and the role of relativistic corrections there. This is what we 
have today under IAU scientific governance over ephemerides centers such as 
the Jet Propulsion Laboratory (JPL), USA, Institute of Applied Astronomy, 
Russia, and some others. 

As for the "anomaly" problem, we are not saying that much better clarity 
than at Clemence's time has been achieved. Some astronomers have contin- 
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ued the anomaly studies basically with the old methodology, which includes, 
on the one hand, a calculation of the Mercury's perihelion advance with "per- 
turbing" forces from all other planets; on the other hand, a comparison of 
theoretical results with empirical data related to Mercury's motion in order 
to fit theoretical ephemerides as close as possible to the Einstein's predicted 
value. 

1.2 PPN and current status of the perihelion advance 
problem 

Work on GR (PPN) extension to the N-body problem was actually started 
by Einstein with coauthors [8] (1949) and was continued during Clemence's 
time and later on, see, for example, [9] (1972), [10] (1973), [llj (1989), [12] 
(1991)), [13J (1993), [H] (2000), [T5] (2003), [T6] (2005), [T7] (2005), [TS] 
(2006), [IS] (2007), [2D] (2008), and references there. 

In the 70s and later on, the N-body problem was theoretically formu- 
lated in the parameterized post Newtonian (PPN) formalisms at the level of 
ephemerides calculations. The GR perihelion advance problem is supposed 
to be included in the formalism. As well known, the GR equations does not 
provide exact solutions for most of practical problems, the N-body system, 
in particular. The PPN approximation idea is to linearize the equations un- 
der weak-field conditions for approximate N-body solutions. Inevitably, the 
relativity essence such as a concept of proper- versus- improper quantities has 
to be sacrificed. 

Our criticism of the PPN methodology is expressed in the question: ap- 
proximation to what? Definitely, it is not the approximation of the exact 
GR solution of N-body problem, because such a solution does not exist. 
This is why the existence of the PPN frame with the corresponding coordi- 
nate system cannot be justified, it is postulated. In the special important 
case, N = 1 plus a test particle, the exact solution does exist; this is the 
Schwarzshild metric, which is valid in the entire range of field strength. In 
PPN methodology, the Schwarzshild metric is replaced with the approximate 
(weak-field) solution containing the PPN parameters /3 and 7. The PPN for- 
malism is intended to account for perturbation of the planetary system due 
to planet interactions, basically, with the use of Le Verrier's method of New- 
tonian mechanics. Nowadays, however, the N-body problem has an exact 
computer-supported solution in the Newtonian model. 
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Table 1: Sources of Mercury perihelion precession 



Amount (arcsec/ century) 


Cause 


5030 


Precession of the equinoxes 


530 


Perturbation by planets 


43 


General Relativity 


5603 


Total 


5600 


Observed 



The incorporation of the PPN formalism into the ephemerides systems 
took place in a process of continual "fitting" of ephemerides systems to ob- 
servation database. Consequently, the PPN parameters were (allegedly) well 
pinned. Overall, an optimistic picture is claimed that astronomical constants 
and planetary ephemerides are adjusted to hundreds of thousands of different 
observations (astronomical, radiometrical and others). Allegedly, any orbit 
can be accurately (in a sense of "fitting" ) displayed in a desirable coordinate 
system and in time directed into past or future for a long enough period. 
We put the word "allegedly" , because there is a difference between "fitting" 
and "evaluation of statistical significance of observational data" , - the issue 
discussed later. 

A practical coincidence of prediction and observation is demonstrated in 
Tabled] (the data in literature usually referred to). Here, we omit originally 
posted precision numbers. In terms of relative accuracies, they are 10 -5 for 
the equinoxes precession (noted by Clemence as the biggest problem), and 
10~ 3 for the predicted GR effect. 

In [T7] (2005) we read, for example: "As the uncertainties in [ephemerides] 
parameters decrease, the domain of possible values of the relativistic param- 
eters narrows, imposing increasingly stringent constraints on the theories of 
gravitation alternative to General Relativity." 

The above and many similar statements should be understood as a claim 
of a huge leap achieved after Clemence 's time in Astronomy with tremendous 
consequences in Physics. There is a firm consensus among astronomy and 
physics communities, mass media as well, that the perihelion advance test 
is the accomplished task. We think, however, that physical reality is more 
complicated, and the above precision numbers for physicists must seem to be 
to a great extent fictitious. 

As a result of our investigation of the problem, we come to the conclusion 



5 



that the long-standing efforts (starting from Le Verrier till present days) to 
confirm the Einstein's prediction of the perihelion advance have not brought 
physical evidence of the confirmation in rigorous terms of physical and sta- 
tistical theories. Further, we present arguments concerning two sides of the 
problem: 

1) Astronomical observations and their "matching" with the Einstein's 
prediction; 

2) Theoretical rigor of the Einstein's prediction. 

2 Calculation versus Observation 
2.1 What and how it can be observed 

Let us see the example of Mercury's ephemerides calculations in the gelio- 
centered frame, Narlikar and Rana [21] (1985), Rana [22] (1987), Figffl A 
classical motion of Mercury's perihelion advance caused by the interaction 
of planets gives the effect 530 arc-sec/per century, see Table [TJ This is a 
fluctuating part of the background, on which the GR predicted effect of rate 
of 43" per century must emerge. 

In the authors' works commented below, a powerful (Burlirsch-Stroer) 
numerical method of differential equation solution is used with the time-step 
about one day chosen to provide a calculational (accumulated in integration) 
error of one part in 10 13 . Such a precision is required to ensure a stability of 
numerical integration over about 300 year period and reverse back. Changing 
the step from 1 day to two days results in a computational error about 0.001 
arc sec. At the same time, initial conditions are known not better than seven 
significant digits. (A position and a velocity of a center of every planet at 
a given instant are meant). This must translate into a cumulative error in 
calculations about 10 arcsec per century. In order to maintain numerical 
stability of calculations, the authors have to control more significant digits. 

As we understand, the computational "instability" in lower order figures 
reflects physical reality of the system of interacting planets (as in Figfl]). 
A small perturbation acquires an ability to propagate into nearby area and 
inflate in time. Indeed, a motion of every planet has many degrees of freedom 
so that an orbit can differ from the elliptic one in however many modes. The 
authors do not conduct calculations in the relativistic (PPN) model because 
of a drastic rise of computer time. This is quite unfortunate because one 
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□ays measured from JD 2444600-5 Number of days following JD 2445400-5 

Figure 1: Position of heliocentric longitude of perihelion of Mercury with time. 
Vertical scale: cumulative effect in arc-seconds. Horizontal scale: days past start- 
ing from December 27, 1980 till April 10, 1984 (left picture [21]), and March 7, 
1983 till February 9, 1988 (right picture [22] (1987)). It is seen that the graph for 
the last 400 days in the left picture and the graph for the first 400 days in the right 
picture overlap in time (as indicated by vertical lines). Total number of Mercury's 
orbital periods is 25. The authors' comments are given in the text. 

would like to see in the PPN formalism the effect of ephemerides fluctuation 
and time instability in the determination of a small regular relativistic effect 
on top. 

The question arises how to determine the difference between "calculation" 
C(T) and "observation" 0(T) for different instants T in view of fluctuation 
of angular motion. In practice, observational data are treated with the use 
of polynomial "averaging" interpolation such as Ti(T) with the term 5600 ■ 
T and To(T) with the term 530 ■ T, so that one has to deal with the problem 
of reduction Ti — > Tq. The authors note that an uncertainty about ±1" 
of geocentric observations translates into the GR effect precision about 10 
%. From the authors' work it follows that the fine structure in calculations 
and observations must be accurately accounted for in order to preserve true 
physical meaning of measured quantities. 

The important part of their works is a criticism of the method of "Ein- 
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stein's effect fitting" used by Clemence and other authors in later times. In 
particular, the results obtained by Clemence, also by Morrison and Ward 
[23] (1975), are discussed. The matter is that the Clemence's results were 
reanalyzed in [23] in view of the fact that precision of data for masses of 
planets and time scalings were significantly improved since Clemence's time. 
This noticeably changed the C(T) while the 0(T) was only slightly cor- 
rected. Nevertheless, the fitting (least square) criterion remained excellent. 
The conclusion of Rana's work [22] (1987) is given below. 

"The low precision of the geocentric angular data having an error of 1" 
are incapable of giving the rate of motion of the perihelion of Mercury to 
better that 3" per century. Hence the determinations apparently good to 
0.3" per century are spurious. Using too low a rejection level only those data 
that happen to agree with the initially assumed value of the rate of motion of 
the perihelion contribute, and one will always obtain a very small correction 
to the initially assumed value. So all the existing observational estimates, 
namely those made by Clemence, Morrison and Ward, and the JPL groups, 
are suspect." 

The author noted that the ephemerides uncertainties must be increased 
substantially. Similar opinion was expressed by Pitjeva [16] (2005). The 
different approach to the problem was made in [24J (1993), when the least 
square criterion was evaluated without paying attention to the Einstein's 
predicted value. After discarding few points, which looked suspicious, the 
author found a drastic drop in the least square number for the advanced 
perihelion value about 15-16" (instead of 43) per century. His comment was: 

"The only conclusion one can draw from the data is thus that they do 
not contribute to a decision as to whether the actual motion of the ascending 
node of the orbit of Mercury exceeds that predicted by the theory" . 

2.1.1 What has been observed 

Now we state that the Mercury's relativistic effect has never been directly 
observed and even not evaluated from circumstantial astronomical evidence. 
The matter is that the GR theory, at least as it given in literature, does not 
provide a clue about distinguishing between the classical drag along with the 
equinoxes precession, on the one hand, and relativistic effect, on the other 
hand. There is no other way but look for an admissible anomaly gap to be 
filled with the predetermined perihelion advance of 43" per century as tight 
as possible, no matter of what kind the effect is. Such "a gap fitting" cannot 
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be termed "the confirmation of the GR prediction" . 

Another issue is a statistical meaning of the gap fitting. There is no single 
publication devoted to the treatment of observations of Mercury perihelion 
advance; the claimed numbers are stated in different works on empirical data 
not treated in rigorous terms of statistical theory. A bad fitting practice and 
the precision concept abuse should be noticed. At the same time, the usage 
of standard "precise" initial conditions in ephemerides calculations makes 
the results stable what creates an illusion of their high-precision, while their 
real precision remains unknown. 

To avoid any terminological ambiguity and fruitless disputes on this im- 
portant issue, we present a brief review of the statistical method that is 
usually used in physical experimental studies and, in our opinion, should be 
used in the perihelion advance investigation as well (see the Appendix A with 
comments). In the following sections, we investigate theoretical rigor of the 
GR predictions. 

3 Classical and GR orbits 

3.1 Classical basic equations 

Before making a comparison of classical (Newtonian) and relativistic (GR/SR) 
formulations of the problem, let us discuss the classical equation of motion 
for a point particle in a spherical symmetric field [25]. The conserved total 
energy eo in the dimensionless form is given by 

eo = i-^ + k 2 + ^4 a) 

r 2 2r z 

where the conserved angular momentum is Iq = rfig = r 2 (d6/dt); r g = GM, 
(3 r = (dr/dt), the speed of light at infinity c = 1. One can replace the equa- 
tion (CQ) with the almost equivalent one having an "relativistic appearance" 
of the squared total energy under the weak field conditions: 

4 = i-& + fi + & ( 2 ) 

Here, the second-order terms are neglected, one of them being Tgl^/r 3 , similar 
to that responsible for perihelion advance in GR. 
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With the use of relationship f3 r = (dr /du)(du/d9)ue, u>e = lo/ r ' 2 with 
u — 1/r, the equation in the standard form is given by 

{du/ddf = (e 2 - I) /ll + 2(r g /l 2 )u - u 2 (3) 

In the traditional form in polar coordinates, the equation has two parameters: 
the eccentricity e > 0, and the so-called semi-latus rectum p (the distance 
from a focus to the ellipse): e 2 = 1 + (ejj — V)r 2 g jl\\ p = l\jr g , 

(du/d9) 2 = (e 2 -l)/p 2 + 2u/p-u 2 (4) 

with the roots: 

u 1 = (l-e)/p, u 2 = (l + e)/p (5) 

leading to 

9 2 - 9 2 = f * du [{u - Ul )(u 2 - u)]- 1/2 (6) 

J Ul 

r{6) =p /[(l + e sin{e-e')\ (7) 
where 9' can be fixed in the initial conditions. 

3.2 Orbit classification: e-criterion versus cr-criterion 

The traditional classification of orbits is based on the e-criterion: 
e > 1, e > 1, hyperbola; 
e = 1, e = 1, parabola; 
e < 1, e < 1, ellipse; 
e = 0, eo = 1 — r g /2r , circle. 

The above classification, though customary in many applications, makes 
a study of the relativistic perihelion problem practically impossible. The two 
parameters p and e are combinations of physical quantities, as they histor- 
ically came from geometrical studies of curves. In Physics and Astronomy, 
one has to deal with a family of orbits specified by one or two physical param- 
eters. In this work, we specify the initial perihelion condition: /3(r ) = /?o, 
at r = r , (9 - 9') = with l = r(3 e = r /3 . 

Further we are going to exploit a non-conventional family of one param- 
eter a = r g /r (3Q for the equation of motion with a radial function being 
= r u(9). In studies of relativistic perihelion advance, a consideration 
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of motion in terms of the azimuth angle 9 was found quite useful. Thus we 
introduce the classical equation in the nonconventional form: 

(di/de) 2 = {l-2o) + 2ai-e (8) 

It is a classical harmonic equation, which can be presented in different forms, 
for example 

d 2 £/d6 2 = <r - £ (9) 
or, by a substitution £ — a = x: 

d 2 x/d6 2 = -x (10) 

The solution to dHJ may be chosen in the form 

e = |H| (ID 

£(0) = r /r(6) = a + (1 - a) cos(#) (12) 

or 

r(6)/r = [a + (1 - a) cos^)]" 1 (13) 

As wanted, the equation is left with only one parameter a, which absorbs r g , 
ro, and (3q. The parameter is intimately related to the 3d Kepler's law. 

Thus, we suggest the following a-classification of orbits illustrated also in 
FigH 

< a < 0.5, hyperbola; 
a = 0.5, parabola; 
0.5 < a < 1, overcircle ellipse; 

= 1, circle; 

1 < a < oo. subcircle ellipse. 

In accordance with the a criterion, eccentric orbits can be of two types. 
The "overcircle" orbits are characterized by initial kinetic energy e& and the 
total energy eo being greater that in the circular motion with r = r , while the 
speed at apohelion is f3 a < (3 . For the "subcircle" orbits it is just opposite: e k 
and eo are less than that in the circular motion, while j3 a > j3 . The subcircle 
motion includes orbits with a — > oo, as f3 — > 0, when the particle's speed 
near the point source becomes however high. This case must be considered 
in the relativistic rather than classical framework. 
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Figure 2: Example of a family of classical orbits. 1. Subcircle ellipse, a— 1.9; 
2. Circle (thick line), a—1; 3. Overcircle ellipse, <r=0.6; 4. Parabola, a— 0.5; 
Hyperbola, a— OA. The gravity center is at the coordinate origin. All orbits 
are produced by launching a point test particle at the perihelion point x—1 
with the initial speed (5q = Jr g /a. 
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In a family of orbits of different radii r (like in a Solar planetary system) 
it is convenient to use the a- classification for the pair of variable parameters, 
ro and 0o, with r$ being the perihelion (minimal) distance from a focal point. 
In this case, the notion of "subcircular" orbit is not relevant. 

The e-criterion cannot, in principle, distinguish between the two types of 
ellipses (overcircle and subcircle) because it suggests variations of physical 
parameters all at once with e. If e varies with p fixed, the family consists 
of orbits with gravity centers sliding along the major axis, what makes a 
physical analysis difficult. Consider, for example, an orbit with a small ec- 
centricity. The equation fll]) has two roots: u — uo(l ± e), or r = r (l =F e). 
Physically this makes a false picture of a point particle being at the same 
time in orbital motions and in a state of radial harmonic oscillation in the 
parabolic well. 

In reality, the particle undergoes an asymmetric periodic motion with 
respect to the gravity center between a perihelion point r = r$ of a lowest 
potential and the apohelion point of a maximal potential. The correct phys- 
ical solution comes from ([H]). The first root is £1 = 1, or r\ = r . The second 
root is £2 = 2er — 1. For a — 1 =p Act, and small Act, it is r 2 = r (l ± Act). 
This gives a subcircle or overcircle ellipses. 

The example of two roots in the ct versus e parametrization of orbits is 
important in view of the idea known from literature about the assessment 
of the GR perihelion advance. It suggests that the orbit slightly different 
from the circular should be analyzed on subject of its angular and radial 
frequencies with respect to the circle frequency. The radial oscillation is 
expected to be the harmonic one in the parabolic potential well. It is thought 
that by finding the difference between angular and radial frequencies one 
could find the resulting perihelion advancement. As shown above, this is 
a wrong approach. We shall return to this issue later in more details. In 
our further analysis of the conditions providing the relativistic perihelion 
advance, the cr-criterion in perturbations of a circular orbit is used. 

3.3 GR Schwarzschild Dynamics of particle and pho- 
ton 

GR Physics of a test particle and a test photon is described by the Einstein's 
field equations [26], [27], or equivalently, the Schwarzschild metric [28], [29J. 
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The latter in polar coordinates is given by: 



ds 2 = dt 2 /r 2 r - r> 2 - r W 



(14) 



where we denote T 2 = 1/(1 — 2r g /r), c = 1, and ds = dr (c = 1). 
The equation is obtained from the Hamilton's variational principle for the 
Schwarzschild field, 1301, I3T1. 



(15) 



5 J ds = 5 J Ldt = 0, L = dr/dt 



with the following integrals of motion 

e = (dt/dr)/T 2 r (16) 

l = r 2 (d6/dT) = r(3 e (17) 

Here, e is treated as the total energy, and 1$ - the angular momentum. They 
are readily follow from the time symmetry and the space isotropy in a manner 
of classical mechanics. The third equation is a definition of the proper time 
( !T4|) usually given by 



{dt/drf/Tl - {dr/dr) 2 Y 2 r 



r 2 {d9/dr) 



(18) 



Here k = 1 for the particle, and k = for the phton. 

The first two expressions are valid for the particle and the photon. The 
equations of motion of the particle and the photon in the Schwarzschild field 
follow from (TT6"j) . f)17p . (lisp . The term which causes the GR effect (further 
called the GR-term) is inclosed into the framed box. The equations will be 
used in our further analysis of GR effects. Below, the orbit equations are 
given 

1. Particle (planet) 



{du/ddf = -(1 - el)/ll + 2{r g /ll)u - u 2 + 



2r g u 3 



(19) 



2. Photon (light ray) 



{du/dOf = el)/ll 



u 2 + 



2r g u 3 



(20) 
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3.4 Introduction of the GR-term in the particle orbit 
equation 

In classical mechanics the "effective" potential energy is 

V(r) = e 2 -(3 2 = l- 2r g /r + l 2 /r 2 (21) 

where (3 2 is the radial part of kinetic energy /3 2 = (3 2 + The extremum 
equation gives a circular orbit with the radius r c 

r c = ll/r g = r g /(3 2 , a = r c /r c ^ 2 = 1 (22) 

The GR perihelion advance is associated with the GR-term (tqIq/t 3 ), 
which is of the second order of smallness {r g ll/r 3 ) ~ i r g/ r ) 2 i n the energy 
balance. At this point, we do not comment the fact that the GR problem 
is formulated in terms of proper time r. It should be also noted that the 
GR theory is supposed to explain laws of gravitational physics exclusively in 
terms of curved spacetime; the concepts of potential and kinetic energy are 
not in the GR arsenal. 

It is habitual to consider the GR-term in the "GR effective potential 
energy" to be compare with the classical case (T2TJ 



V ef (r) = 2r g /r + l 2 /r 2 - 2r g l 2 /r s (23) 



In the angular part (3$, the quantity Iq = r (3 must be the conserved orbital 
momentum. 

The equation fl23l can be rewritten 



Vef /p 2 = l/ffi-2*Z + e- 2(r g /r )e (24) 



Finally, it can take the form of equation of orbital motion, by analogy with 



(dt/d9) 2 = (l-2a) + 2at-e+ 2(r 9 /r )£ 3 (25) 



The above GR equations have common extrema and can be used for deter- 
mination of orbital motion parameters. The extrema are readily found by 
equating the derivative of right side of equation to zero. 

Let us see potential energy curves with and without the GR-term for a 
strong field condition in order to understand the role of the GR-term. The 
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graphs (I2"3"j) for the classical potential energy (upper line, no GR-term) and its 
GR version (below, with the GR-term) are shown in FigJHJ Both curves are 
plotted for Zq = 20 in r 2 g units. Radii are given in units of r g . Consequently, 
the classical radius of circular motion is r c = 20 and the corresponding speed 
/3q = 0.22. For the GR curve the radius of the circular motion is about 
r = 16.8. The "peak" at the radius r\ = 3 is caused by the GR-term. 

Here is "an inner" radius about r* = 2.5, which, in our opinion, is a 
fictitious solution. One can consider it a solution for a confined particle in a 
separate circular orbit slightly above the so-called Schwarzschild sphere with 
a speed of motion close to cq. If to take r g = 1.5 • 10 3 m (the gravitational 
radius of the Sun), the orbital radius will be R = 2.3 • 10 4 m. Compare it 
with Mercury parameters: semi-latus rectum p = 5.8 • 10 10 m, (3 = 1.3 • 10~ 4 . 
Looking at the picture, one has to bear in mind the GR assertion that the 
GR-term drives the perihelion precession in the direction of motion. 

3.5 Classical versus GR treatment of the GR-term 

The GR-term cannot arise in the Newtonian orbit equation except for usage 
of sticks and ropes. Having this said, one can treat the equations with the 
GR-term classically (consistently with the criterion a), as in (1251) . For r$ = 
r c (l — 3r g /r Q ), we have a = r g /r c (l — 3r g /r c )f3 2 > 1 that shows a subcircle 
ellipse. It can be transformed into a circle by variation of 0: with 0i = 
0o(l + 3r s /r ), we have a = 1. 

In our approach, the cubic algebraic equation has always a unit root 
£_ = 1. One can find the circular solution or the elliptic ones along with 
the third root £* = r /r* to the precision of any higher degree and reduce 
the equation to the quadratic form. The main solution is just an ordinary 
classical one. One can vary physical parameters either keeping Iq constant or 
redefining it. The main thing to remember is that the classical circular orbit 
always abide by the inertial rotation law a = 1. 

Suppose, from physical considerations, it is known that the circular orbit 
is characterized by parameters, which make the criterion a ^ 1. This would 
be a sign to treat the equation of motion alternatively, say, in terms of SR or 
GR. Such an example is known in SR Kinematics, - the Thomas precession, 
discussed later. 

Now let us consider the approach to the GR perihelion advance problem 
based on the GR concept of effective potential. The latter can be "con- 
structed" from the Schwarzschild metric, as discussed in the famous mono- 
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Potential energy 




Figure 3: GR "effective" potential energy (thick line) as a function of radius 
in r g units. A circular motion radius is r . The upper curve is the Newtonian 
potential (the GR term removed) with the corresponding radius of circular 
motion r c . In both cases, Zg = 20 r 2 g . Notice the shift of two graphs due to 
the GR-term. For comments on r* and other details, see the text. 



graph "Gravitation" by Misner with co-authors [10] (1973). They suggested 
that the Einstein's perihelion advance of 3r g /r (rad per 2ir) could be con- 
sidered in an approach different from the Einstein's one. The starting point 
is to consider an orbit with a slightly different energy than in a pure circular 
motion. The authors gave some hints to the solution but suggested that stu- 
dents would do the work. Such a work was actually done by R. Wald (1984) 
[32] . and it was widely popularized later, for example, [29], (2000). The as- 
sumption (mentioned above in connection with the e-criterion of orbits) is 
made that the value of perihelion advance could be calculated as a difference 
between the angular frequency uig and the frequency of radial oscillation uj t . 
Due to the GR-term impact, the angular frequency increases by the factor of 
(1+ 1.5 r g /r ). The radial frequency must be found from the second deriva- 
tive of the potential as an approximation of the law of harmonic oscillation 
in a parabolic potential well. One can visualize the picture of harmonic os- 
cillation by examining the solution uq ± e from the e classification of orbits. 
The radial frequency must decrease by the factor of (1+ 1.5 r 9 /ro). A to- 
tal relative change of a classical angular frequency by a factor (1+ 3 r g /ro) 
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appears to be in agreement with the Einstein's prediction. 

There are objections. First, a comparison of classical and relativistic cases 
needs relativistic- versus-classical arguments and criteria. Second, in order to 
make the perihelion advance in the direction of motion, the angular frequency 
must decrease (not increase). We shall discuss this issue later in more details. 
Third, the radial and angular oscillations cannot be manipulated in arbitrary 
perturbations; they are engaged, as seen from relations (3 r ~ uj r ~ dr/d6 ~ 
ug. Fourth, the picture of parabolic well appears in the e orbit classification, 
which is, though geometrically correct, physically misleading and wrong in 
this discussed earlier. 

We think that work with the relativistic equations ( I23p . (I24p . ( 12 5 p in 
Newtonian environment is, in principle, a wrong idea. Einstein's original 
approach was technically and methodologically different. He worked with 
the equations embedded into the GR framework with conceptual attributes 
such as symmetry and conservation laws, particle and photon parallels, fi- 
nally, proper versus improper time scales. Whether he took advantage of 
it is another matter. Anyway, one can raise the question why to search for 
alternatives when the solution of the perihelion advance problem was given 
by Einstein in 1915 and never seriously doubted. Or was it? 

4 The GR angular advance, and the SR Thomas 
precession 

4.1 The SR Thomas precession 

The mentioned earlier Clemence's notice that the GR perihelion advance is a 
rotation of orbital plane deserves a full attention. Most likely, Clemence was 
aware of the SR Kinematics effect, - the Thomas precession [33] (1926), which 
is only one sixth the GR effect but resembles it in all other respects. Can the 
Thomas precession contribute to the GR effect? To clarify this question and 
find out the ways of distinguishing between classical and relativistic motion, 
we need to discuss the phenomenon. 

In a rotating frame, centripetal acceleration is caused by a stress, which 
affects clock rates. This suggests that the Thomas precession has to be 
considered in terms of SR Dynamics. Nevertheless, it is often treated in 
terms of SR Kinematic as an approximation tolerable to the extent when 
dynamical corrections to the effect are neglected. Still one has to bear in 
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mind that in SR Kinematics of uniform linear motion there are no preferred 
coordinate system, while a system uniformly rotating by inertia is a non- 
inertial, hence, preferred one. 

Consider a uniformly rotating rigid frame with a coordinate system S' 
attached to the frame. The system rotates about the origin O' with the 
angular speed u . The rotation is seen from an inertial coordinate system S 
with the center O coinciding with the center O' both being at at rest with 
respect to the stars far away. Notice, radial scales in S and S' are similar. 
A point r' = ro is in a circular motion in S; it moves along the circular arc 
s with the speed (3q = r u. The arc path interval is As(t) = PoAt = r uAt, 
where t is a coordinate time related to the proper time of rotating observer 
t' = t/7. The Lorentz factor is 70 = (1 — /3q) _1 ^ 2 , hence, s/s' = 70 . This is 
the SR Thomas precession phenomenon. 

4.2 The panetary version of the Thomas precession 

Let a hoop rotate about the gravity center at O coinciding with O' so that 
the a-criterion a' = r g /r f3Q — 1 in S' and, similarly a = Tgjr^pj^^ = 1 in 
S is satisfied. In both systems, the hoop is subjected to no internal stress. 
For a > 1 the stress would arise due to hoop stretching, for a < 1 for 
compressing. Now, cut the hoop into small pieces, and we have the version 
of Thomas precession of a planet in a circular motion. The orbit will acquire 
an elliptic, nearly circular, form if a small instant perturbation 5r is given. 
Hence, the arc advance can be observed in S in the form of precession of 
the line connecting perihelion and apohelion points. In both systems, the 
following definition of the half-circle period applies: the time passed along 
the arc nr from the event "perihelion" to the event "apohelion" in a limit 
5r ->■ 0. 

A planet or satellite precession can be observed in a pure circular motion. 
If the clock in S' "flashes" every period, the flash point is seen precessing 
in S in the direction of motion. The comoving observer in S' cannot see 
the precession unless she periodically snaps stars far away. Mathematically, 
the orbit precession is described by a periodic function with the argument 
ut/j ps vut, where v = (1 — /5q/2), (3q = r g /r , in accordance with the 
criterion of a — a' — 1. The reduction of the angular speed makes a margin 
to extend the period in S' to its value in S. We think that the SR Thomas 
precession must take place and can be observed in a planetary motion. 

The Thomas precession is loosely termed "the orbital plane rotation in 
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the direction of motion". However, such a terminology belongs to classical 
physics language, which in no way recognizes the Minkowski space concepts 
including Lorentz transforms S' # S and rescaling of physical units. 

4.3 About relativistic terms 

Many disputes about relativity (not to speak about GR) are counter-productive 
just because of a loose terminology, an abuse the SR concepts and definitions, 
finally, neglect of operational meaning of the Lorentz transformations. See 
some examples below. 

Relativity. One should bear in mind that the SR concept of "relativity" 
of motion cannot be "refuted" because it stands for the classical (Galilean) 
relativity, on which the SR is based (the postulate 1). The SR Theory cannot 
be refuted because it is based on the postulate 1 and 2, the second postu- 
late being a constancy of speed of light in all inertial coordinate systems. 
Both postulates are firmly proved in experiments. Why the GR is " General 
Relativity" is actually not clear, see [3Tj . 

Event. A non-trivial "event" includes two events constituting the 4- 
coordinate and 4-momentum vectors transformed by the Lorentz operator. 
A vector must have the origin and the tip, the two events, which are: 

1. synchronization of clocks in S' and S at t' = t = 0, x' = x = 0; 

2. light signal exchange between S' and S observers. The concrete pro- 
cedure depends on an objective of observation (type of intended event). 

Coordinate time. There is no such a thing as an observer's wristwatch 
recording the coordinate time t: the time is measured with the use of records 
of S"-clock at x' and S-clock at x. In "the pure Lorentz boost", y and z 
coordinates are not involved. The time dilation and the length contraction 
take place only in the direction of motion. A motional (angular or linear) 
advance is an "inverted" length contraction. 

Clock paradox. The problem with the paradox typically lurks in the ill- 
posed formulation and a potential abuse of the coordinate time concept. The 
matter is that a rest observer can work by herself "observing" events from 
records of proper time intervals. Alternatively, she can be in couple with 
the inertial observer to collect data for measuring the coordinate time. The 
quantity to be compared is the Lorentz invariant one such as the proper 
time interval, or the phase of periodic events that is, a number of clock ticks 
N detected by a wristwatch of each observer. A different approach to the 
paradox formulation is a collection of periodic light pulses in the Doppler 
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operational mode [31]. The clock paradox can be particularly recognized in 
the Thomas precession problem. 

Further we continue to analyze the perihelion advance problem. So far, 
we failed in our efforts to realize the physical meaning of the GR-term at a 
deeper than algebraic equation level. We are looking with hope to finding 
answers in original works by Einstein and prominent GR experts. After that, 
we shall return to the question whether the Thomas precession can be seen 
in the sky as a part or the whole relativistic gravitational effect. 



5 The GR-term with Einstein 
5.1 The year 1915 

Einstein starts with the equation ( I19p . [26] (1915), which is derived in the 
scale of the proper time r. Still it is not clear why r but not t. He treats the 
solution in terms of e-representation; hence, the equation is given by 



(du/dOy = -(1 - e l )lv + (2/p)u - vT + 2r g u A (26) 



where u — 1/r, e and p are classical ellipse parameters. An elliptic solution 
has the form u 1)2 = u (l ± e). The equation has three roots, %, it 2 , u 3 , 
where r\ and r 2 are perihelion and apohelion points while r 3 has no classical 
analogy; its origin is due purely to the GR term, see the root r* in Fig. |3J 

The Einstein's plan to solve the perihelion advance problem was straight- 
forward : to integrate the arc path from perihelion to apohelion to see how 
the result differs from tt. Many years later, the plan was reproduced by Moller 
[35] , (1972) with somehow different physical concepts. He came to the same 
integral, as Einstein did, but in the scale of coordinate time t. Anyway, the 
solution is 

/"^2 r 1—1/2 

9=1 du (e 2 - l)/p 2 + {2/p)u - 2r g u 3 (27) 

It should be noted that the Moller's derivation is presented in more details. A 
reader is recommended to read both works in parallel (Einstein's in German 
and Moller's in Englsh). Notice a missed factor 1/2 in (16) in [26]. 

Sadly, the Einstein's and then Moller's plan was not accompanied with 
discussions of the role of the GR framework and the corresponding initial 
conditions. As shown above, without comparing circular motion conditions 
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in classical and relativistic variants (say, by a criterion), the problem be- 
comes ill-posed, and a conclusive result cannot be achieved. The same re- 
marks can be addressed to the Synge's work [36J. Our critical argument on 
the methodology of the GR effect based on understanding of the fact that 
the relativistic angular advance (in both SR and GR methodologies) has no 
analogy in classical physics. The relativistic criterion of observability and 
treatment of the effect should be described in clear GR terms. Let us see 
how the plan proceeds. 

To avoid complications with the exact solution (I27p . Einstein considers, 
as an approximation, the following integral 

1/2 

ill 



6=1 dw[-2r ff /r («-Wi)(w-w 2 )(«-W3)] 1 (28) 



where u\ and u 2 are the exact roots for ([19]) . Then he decides to modify the 
polynomial of third degree in order to use the exact solution of the classical 
equation ui and u 2 fl5]) assuming that the equalities u\ = ui, u 2 = u 2 are valid 
to the precision of smallness of the GR-term. In other words, Einstein (and 
later Moller) makes the assumption that the main roots u\ and u 2 are not 
affected by the presence of the third root u* = 1/r* so that u* » {ui + u 2 ) 
to justify the assumption {u* + u\ + u 2 ) = l/(2r g /r ) for however small e. 
This makes the third root a known function of the main roots u* = f{u\, u 2 ) 
so that the multiplier v — (1 — 3r 9 /p) appears in the argument of periodic 
function 



(l/p 2 )(du/dey = (e 2 -l)/p z + (2/p)u-u 2 + 2r g u 3 (29) 



the solution of which describes the elliptic motion with the argument (vO): 

u = l/p+(e/p)cos{v6) (30) 

We argue, however, that the Einstein's assumption is physically and log- 
ically inconsistent: if the third root u* (which is closely related to the GR 
term) does not affect the main roots ui and u 2 , then it cannot be presented 
by a certain function of them. In fact, we know that the shift of main roots 
is about 3r g what is comparable with the GR-term and the 3d root. 

We state that the Einstein's wrong assumption and rude approximations 
lead to the false prediction of the perihelion advance 3r g /rQ. It can be ver- 
ified by a direct comparison of the original algebraic expression with the 
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approximate one obtained by Einstein [25] (1915) 

- 2r g [(it - Ul ){u - u 2 )]- 1/2 = (1 + 2r g /p){\ + r g u)) [{u - u x ){u - u 2 )]- 1/2 

(31) 

Squaring and inverting both sides of (13"T1) leads to the contradicting ex- 
pression 

= -2r g /p (32) 

Our argument that the Einstein's prediction is not valid can be also verified 
by a direct substitution of (1301) into (1291) . Having this done, one finds an 
exact cancellation of all terms on both sides of the equation but the GR 
term on the right: 

ee I 3(r g /r )u 2 I (33) 



Thus, the claimed formula for the perihelion advance 3r g /ro is a result of 
inappropriate mathematical assumption in the equation solution, therefore, 
the prediction is not valid. 

5.2 The year 1942 

Probably, Einstein did not feel a satisfaction with his account for the GR 
perihelion advance. In the following years he never returned to the above 
derivation in his publications and never suggested some modification of it. 
The paper [26j (1915) was not translated into English, and still it remains 
hardly accessible. About a quarter of century later, the monograph on Rel- 
ativity Theory by Bergman [27] (1942) was published where Einstein in the 
foreword acknowledged his personal role in its composing. There, for the first 
time another derivation of the effect is presented. It starts with the second 
order differential equation 



du 2 /d6 2 = l/p-u + 3r g u 2 (34) 



where the GR term carries the coefficient responsible for the main root shift 
and the predicted effect. Instead of integration of the equation, the author 
suggests that the Fourier's expansion of the solution be analyzed in view of 
a small parameter A = 3r g /ro (in his denotations): 



u 



ao + ai cos (p9) + A v 2 /3 u cos (vpQ) (35) 
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starting with the first approximation 

u = a(l — e cos)6> (36) 

After making some more (not clear) approximations, the author comes to the 
solution obtained before ( 1301) . which we have already proved to be wrong. 

5.3 Fock's work, 1955 

About 13 years after the Bergmann's work, the monograph on "The Theory 
of Space, Time, and Gravitation" by Russian physicist V. Fock was published 
with a critical attitude towards the GR theory [31], (1955). There Fock sug- 
gests the concept of harmonic coordinate system, the idea of which reflects 
d'Alamberian wave conditions at infinity. Thereafter, the Schwarzschild met- 
ric undergoes some modification with the appearance of terms of higher (up to 
fourth) order in the equation (1191) . While considering the Mercury perihelion 
advance, Fock follows the Einstein's methodology of relativistic reduction of 
frequency in periodic motion u = (1/p) (1 — e cos (u9) with v — (1 — 3r g / r ) to 
provide the perihelion advance. His analysis of the problem in the harmonic 
coordinate system shows that a combination of quantities in the right side of 
the equation can match the multiplier v needed for the relativistic rotation of 
the orbital plane. However, the terms of 3d and 4d order must be neglected 
(the Einstein's GR-term included). It means that the GR term plays no role 
in the effect formation. However, the introduction of "harmonic conditions" 
can make a difference. 

5.4 The GR problem of the perihelion advance: sum- 
ming up 

We state that the GR prediction is in a serious doubt. First of all, the problem 
is ill-posed: its formulation does not allow an observer to differ between a 
small GR effect from similar classical effects. Conditions, under which the GR 
effect acquires relativistic properties consistently with the GR framework, are 
not specified (for example, the role of r versus t scaling). There is no answer 
to the question whether the kinematical effect of the Thomas precession 
contributes to the GR effect in question. Finally, technical problems in the 
derivation of the GR effect should be noted. Namely, the constraints imposed 
by the e-classification of the orbit make a physical analysis of the problem 
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ambiguous. The actual derivation of the solution is flawed. Overall, the 
claimed prediction is proved physically invalid. 

Our expectations to clarify the problem formulation from the original 
Einstein's work and works by reputed GR experts over decades after the 
Einstein's work did not come true. In literature, authors present derivations 
of the effect with a variety of different techniques, but always with the same 
result. At the same time, they do not attempt to analyze factors of GR 
spacetime curvature, which would create the effect. No one wanted to check 
the solution by numerical computations or by its substitution to the original 
equation. The hard questions remain if the effect can survive without the 
GR-term (the Fock's work), why the effect in most of the works is derived 
in the scale of proper time r rather than the coordinate time t (the Moller's 
work), finally, what role, if any, plays the GR-term in the GR testable effects 
with light (to examine immediately). 

The above statements, if true, are too serious and raise questions about 
validity of the GR physical foundations. Though it is out of the scope of 
the present work, some issues of GR foundations will be briefly discussed in 
the remaining part of the work. In particular, we want to investigate a role 
of the GR framework in providing a unified approach to problems involving 
particles and photons. 



6 GR tests with light 

6.1 The GR-term as the cause of the bending of light 

Previously, it was shown how the GR equations of motion of the particle and 
the photon are derived from the Hamilton's principle of variation applied 
to the Schwarzschild metric. The two conserved quantities are found: e, 
associated with the energy (flBT) . and l - the angular momentum (fTTj) . The 
remarkable (and strange) fact is that those quantities are the same for the 
particle and the photon 

e = (dt/dr)/T 2 r , l = r 2 {d6/dT) 

This is strange because the angular momentum is given per unit mass, and 
the question arises, what mass the photon has. The equation (1201) describing 
the bending of light in the Schwarzschild field is 



(du/dey = e 2 /r -u 2 + 



2r g u 3 
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Now a reader should be especially attentive. Unlike in the case of particle, 
there is no linear term there related to the gravitational potential. Hence, the 
GR-term solely is supposed to be the cause of the bending of light. Indeed, 
putting the term to zero results in the solution describing a straight line 
trajectory of the photon similarly to the classical photon motion in "empty" 
Euclidean space. Further we shall see grave consequences of the lost linear 
term. 

6.2 Einstein's evaluation of the light bending, years 
1911, 1916, 1947 

The Einstein's first assessment of the bending of light grazing the Sun is 
given in [37], (1911), A9 = 2r g /R, where R is the Sun's radius. Though 
Einstein refers to the Huyghen's principle he actually uses the idea of light 
refraction with the index n r = 1/(3 ((the Snell's law)) to integrate the angle 
along the light path. This means that the speed of light depends on the 
potential /3(r) = (3q(1 — r g /r). The result was widely criticized at that time 
because it was exactly the same as obtained in the Newtonian model years 
before. 

In the second paper devoted to the GR foundations [38], (1916) Einstein 
corrected the calculation. Again he refers to the Huyghen principle but 
makes calculations of the same integral with the redefined index of refraction 
n r w (—gu/g22)^ 2 ~ (1 + 2r g /r). This should double the effect making the 
full angle of deflection (between two asymptotes) 4r g /R, about 1.7". 

As it seen, both earlier and this evaluation has nothing to do with the 
above discussed GR framework. Einstein, when working on the perihelion 
advance problem in 1915, had to be aware of the space-time symmetries, from 
which the equations of motion ( TT8|) and the trajectory for both the particle 
and the photon were deduced f)19p . (I26p . However, his derivation of the light 
deflection in 1916 was made in the ad hoc approach, outside the above GR 
framework. The effect about 1.7" of deflection was confirmed in a series of 
observations of total eclipses of the Sun's, but this confirmation was made 
against the logic of GR physical foundations rather than in accordance with 
that. This circumstance remained largely unnoticed (or thought excusable?) 
in the GR history. 

Undoubtedly, Einstein realized the methodological inconsistency of his 
predictions and the fact that the derivation of the bending of light had to be 
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reexamined and formulated in the unified GR framework given by (I16p . (|17p . 
(ITS]) . dED, (1201) . About fifteen years passed before it was made. In 1942, 
Bergmann (with Einstein) [27] abandoned the idea of refracting deflection 
and finally turned to the equation (1201) in order to evaluate the bending of 
light 



{du/def = el/ll - u 2 + 2r g u 3 (37) 



The approach used this time is consistent with that for the perihelion advance 
problem. But now (recall our warning) one encounters a severe problem 
of absence of the potential term in the photon trajectory equation. The 
only term left to cause the effect is the GR-term. The latter, as we found, 
plays a controversial role in the perihelion advance problem. In particular, 
the result of Fock's work, [31] (1955) is remarkable: the perihelion advance 
is explained due to the Schwarzschild metric adjustment to the harmonic 
coordinate system, while the GR-term is ignored. As concerns the role of 
potential energy term in the bending of light, see also the textbook [TO] 
(1970): 

"Relativistic effects [in GR] on light and radio-wave propagation are gov- 
erned entirely by the Newton potential U and the PPN parameter 7" . 

The matter is that the authors [JO] were not inclined to investigate the GR 
framework; they evaluated all the GR classical effects in the PPM formalism. 
We do not know how the potential term was reborn there. 

In 1942, nothing could be done but take the same plan as in the perihelion 
advance problem: to find corrections to the exact unperturbed solution (a 
straight line for r g = 0) due to the "perturbation effect" caused by the small 
GR-term. Interpolations are used based on the assumption that the exact 
roots of ( 157]) are not influenced by the small term. 

The author [27] arrives at a "correct" result AO = 4r g /R= 1.7" with the 
use of procedure, which is hard to follow. It should be noted that the solution 
was "reproduced" later, for example Moller [55] (1972). In [32], also [2JJ], the 
GR concept of the effective potential V e f is used to get the same equation. 
The PPN formalism was used in several works, [10], [9], and other works. 
They are not discussed here for the reasons explained earlier. 

In order to see the inherent inconsistency of the Bergmann (with Ein- 
stein) solution of the light bending problem, one has to unfold the algebraic 
procedure to the point of checking the result by substituting it back into the 
original equation (as we did in the case of Mercury's perihelion advance). 
We found that the claimed solution is not only technically but fundamen- 
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tally wrong because of the loss of the potential energy term. We need now 
to discuss the Fock's work [21] (1955), which brings a lot of clarity in this 
problem. 



6.3 Fock's work, 1955 

The Schwarzschild metric in Fock's harmonic coordinates (21] is given by 

* 2 = (^) * 2 - (^) *- 2 - ( r + ^ < 38 > 

The photon equation follows from (J3"8"j) after setting up the boundary condi- 
tion and neglecting terms of orders higher than ~ u 2 

{du/dOf = 1/R 2 -u 2 + Ar g u/R 2 (39) 

Compare it with (JSTj) (Einstein 1947, Moller 1972, Wald 1984) 

(du/dB) 2 = l/R 2 -u 2 + 



2r g u 3 



Here R is the impact parameter (the Sun's radius in eclipse observations). 
In the Fock's equation, this is a liner (potential energy) term that makes the 
light bending effect. The solution is 

Ru = 2r g /R + cos Q (40) 

It describes a hyperbolic trajectory. The full bending angle (between two 
asymptotes) is 4r 9 /i?=1.7". 

Fock's comments follows. The effect can be treated in terms of re- 

fraction in the gravitational field when the latter is considered the optical ac- 
tive medium with a continuous change of the index of refraction n r = l//3(r). 
This suggests that the speed of light decreases with the potential depth. If 
the effect is treated in terms of gravitating photon in Newtonian Physics (the 
speed dependence on the potential only), then the effect would be half the full 
one (the value obtained by Einstein in 1911). The photon motion equation 
derived from the Schwarzschild metric in the harmonic coordinate system 
gives the correct effect due to the potential energy term leading to n r = T 2 
that is, /3(r) = 1/r 2 = (1 — 2r g /r). However, usage of the Schwarzschild met- 
ric in "the standard coordinate system" (Tl4"|) results in the motion equation 
fl2"U|) with the potential energy term lost. 
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It looks like the Fock's harmonic conditions make a difference. In connec- 
tion with the choice of coordinate conditions in GR, it is usually emphasized 
the fact that the physical world does not depend on our choice of coordi- 
nate system (the general covariance). This would be a simplification of the 
problem. In GR, not just a coordinate system should be fixed. In order 
to uniquely solve the Einstein's field equations, four differential equations 
are provided that the metric tensor must satisfy. Therefore, the coordinate 
conditions have to impose physically motivated constraints on solutions. 

The Fock's harmonic coordinate system is not original; it is also known 
from the earlier suggested de Donder's gauge: the 4-coordinates, which are 
Lorentz invariant solutions to the d'Alamber's equation. In [3T] , motivations 
for choosing the harmonic conditions are thoroughly discussed, in particular, 
in connection with the conditions at infinity. 

6.4 The Light and the Particle in GR 

6.4.1 Brief conceptual review of the GR tests with light 

In the Schwarzschild field, there is a symmetry between the photon and 
the particle. Expressions for their total energy and angular momentum are 
identical. The only difference is in that one has to put dr = for the photon. 
The photon can be treated as the particle in a zero mass limit still not loosing 
a property of attraction by the gravitational source (because the energy is a 
source of gravitational field in GR). In particular, the photon can be launched 
into an orbit in a particle manner. 

However, the GR photon concept dr — does not define the photon wave 
properties such as frequency and wavelength. A connection of those proper- 
ties with a photon source is not recognized either. Moreover, the de Broglie 
wave concept in the modern Special Relativity theory is not respected be- 
cause the Special Relativity theory is declared incompatible with the gravity 
phenomenon, paradoxically, for an alleged failure to describe physical prop- 
erties of the light, or photon, [TO]. At the same time, we just witnessed the 
GR failure to explain observations of the the light bending effect, where wave 
properties are most likely involved. What about the red-shift effect? 
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6.4.2 The red shift 



Consider the photon emitted from the Sun's surface and observed on Earth 
(the GR red-shift test). The Einstein's explanation of the observation [38] 
(1916) is based on the comparison of frequencies of the standard clock and 
the photon. Let the standard clock have the frequency fi(cl), say, on the 
Sun's surface r± that is /i(cZ) = f (l — r g /ri). Then he notes that the photon 
of the same frequency fi{ph) travels from r± to r2, (say, the Earth) and is 
detected there with the same frequency fi(ph) = f2{ph), while the clock at r2 
runs faster, namely, f2(cl) = fo(cl)(l — rg/r 2 , fo is the clock rate at infinity. 
However, he could give another scenario of the event development with the 
same assurance, as next. A standard photon source is not influenced by 
the gravitational field. This is the frequency of the photon in flight, which 
depends on the potential. This would also explain why the Earth's observer 
detects red-shift in a star light. 

The GR expertise "explanatory" comments on this issue [39] , [9] do not 
and cannot bring clarity because gravitational properties of the photons and 
the photon emitters are not defined in GR. A relationship between the pho- 
ton's gravitational characterisics such as the speed f3(r), the frequency f{r), 
and the total energy e(r) = (dt/dr)(l — 2r g /r) should be brought into the 
GR framework. 

Let us sum up the bending light and the red-shift issues. In the Bergmann- 
Einstein's explanation [27\ (1942) of the light bending effect, the photon is 
treated as a particle having no wave properties, consistently with the GR 
framework ffT4|) . ffTBT) . f|T7|) . ffl~9l) . fl20|) . The Focks' result (under harmonic 
coordinate conditions) also was obtained without any additional postulates. 
The qustion arises why the predictions are different. 

The author [39J seems to be aware of the controversy and suggests to mod- 
ify the GR foundations, namely, to add a postulate connecting gravitational 
properties of a photon source and the photon. Specifically, the frequency of 
the emitter of the standard photon should be a function of the gravitational 
potential f(r) = f (l — r g /r), while the frequency of the photon in flight does 
not change. 

6.4.3 The particle in free radial fall 

As previously shown, the prediction of GR perihelion advance caused by the 
GR-term is not valid. We state that the particle orbit equation, though being 
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consistent with the GR framework, is physically wrong. There is another GR 
prediction rarely discussed in literature. It concerns a particle radial motion 
in the gravitational field. The GR theory predicts that the particle approach- 
ing the center in free fall decelerates, consistently with the GR framework. 
This effect is, in principle, testable, for example, in observations of cosmic 
high-energy particle coming onto Earth. 

Specifically, [10J, also [10], the speed of a particle in the radial fall is 



Here, 7 = Eq/itlq > 1 is the initial total energy at infinity. A free fall from 
rest corresponds to 70 = 1. The formula is given in a coordinate system 
of the observer at infinity at rest with respect to the gravitational center. 
Therefore, her wristwatch time is the coordinate time t. 

The formula shows [40J that the particle sent from infinity to the gravity 
center begins to accelerate, then at some point starts decelerating and even- 
tually stops at rd = 2r g . The higher initial kinetic energy, the farther the 



point of r d from the center. For 7 > y 3/2, the particle will never accelerate 
in a gravitational field. The gravitational force exerted on the particle be- 
comes repulsive in the entire space. We think, however, that there is no real 
physics behind this prediction |41j . 

One should notice that in GR problems the quantity e plays a special role 
in the formulation of initial conditions. Next, a controversy with GR concept 
of energy e = (dt / dr) /Y 2 is discussed. 

6.5 The GR energy and the initial value 

Recall the GR scheme: from the Schwarzschild metric (fl4"|) 



P(r) = (1 - 2r»[l - (1 - 2r g /r)/ 7o 2 ] 



(41) 




dr 2 = dt 2 /T 2 r - T 2 r dr 2 - r 2 d9 2 



with 




{du/dd) 2 = -(1 - e 2 )/l 2 + 2(r g /l 2 )u - u 2 + 2r g u 3 
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Typically, e is absorbed in the orbit equation without need to specify the 
quantity (dt/dr). But often one needs to express the initial values Iq and 
e = [{dt/dr){\ — 2r g /r)] in terms of r , #0 and the squared speed (3q = 
((3 2 + /3|)o in an arbitrary geometrcal configuration. In GR text-books, for 
example [29] and elsewhere, there are instructions on how to do this. It is 
suggested to introduce "the shell observer" who accounts for the proper time 
r s h with respect to the variable proper time r, also, keeps connections with 
the coordinate time t, as next 

(dt/dr) = (dt/dr sh )(dr sh /dr) (42) 

where 

(dt/dr sh ) = (1 - 2r g /r)- 1 / 2 ; (dr sh /dr) = (1 - (3 2 - /3 2 )" 1 / 2 (43) 

Hence 

(dt/dr) 2 = (1 - 2r,/r)/(l - (3 2 - (3 2 ) (44) 

where (3 2 = l\/r 2 . 

As a by-product of the above "instruction" , we have a new expression 



From 



(1 - 2r»/(l - d 2 , - ll/r 2 ) = (1 - 2r,/r)/(l - & 2 - / 2 /r 2 ) (45) 



/3 r 2 = (6 2 - l)/e 2 + 2r,/6 2 r - / 2 /r 2 + □ (46) 



But where is the GR-term? It is gone. Recall, the GR-term appeared in the 
equation of motion from the Schwarzschild metric, and "the instruction" is 
supposed to be consistent with the GR framework. But it is not. If "the 
instruction" has a physical meaning and removes the GR-term from the box, 
it explains why Einstein and other authors suffered so much with the GR 
term. It looks that the GR field equations and the following from them 
Schwarzshild metric have a serious initial value problem. 



7 Different Theory? 

We propose another (relativistic gravitational) theory In the referred 
work, the principles of relativistic dynamics the theory based on are given 
and illustrated in the example of free radial fall in the spherical symmetric 
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field. The orbital motion problem is described in the present work, as we 
planned. 

In brief, the principles £1X6 5 £IS follows. A new concept of the relativistic 
proper mass m(r) depending on field strength is introduced. From the La- 
grangian problem formulation, it follows m(r) = mo/^ r where m(r) — > mo 
as r — > oo, with 7 r = exp(r g /r). The revision of the proper mass concept 
is motivated by several reasons, one of them, a necessity to introduce the 4- 
momentum vector P M in the form complementary to the 4-coordinate vector 
X^. The temporal component in X^ is the proper time depending on the 
gravitational potential r = r(r ff /r). Therefore, the temporal component m 
in P M should be m = m(r g /r. This explains the gravitational time dilation. 

Thus, the gravitational dynamics is formulated in the Minkowsli space in 
the presence of gravitational sources. In polar coordinates, the 4-coordinate 
interval and the 4-momentum vectors are dX^{r) = 7<ir(r) (1, fl r , flg) and 
P M (r) = 7m(r) (1, (3 r , (3e) , where 3-velocity components and the Lorentz 
factor are functions of r and 9, Cq = 1. The Minkowski 4-force = dP^/dr 
acts on the test particle, and it naturally has the tangential component (with 
respect to the world-line s) and the orthogonal one, while s is a function of 
4-position. 

There are two conservation laws, - for total energy eo, and the angular 
momentum L given below for initial conditions r(r) = r , 9 = 0, (3 r = 0, 

fa = /?o 

The total energy and the angular momentum are 

e = 7o7r = 77r (47) 

L = 7o7ro r oA) = llrrfie (48) 
Instead of (14"8"|) . it is convenient to use a conserved quantity l = e /L Q : 

k = r(3 e (49) 

Here, a squared inverted Lorentz factor is I/7 2 = (1 — $1 — /3|), (3 r = 
dr/dt, (3g = r(d9/dt), and we are going to use, as usual, the formula f3 r = 
(dr/d9)(d9/dt), and (3g = l^/r 2 . After introducing a variable £ = r /r, we 
arrive to the exact relativistic equation of orbital motion of confined particle. 
The equation has the Newtonian limit, and it is valid for a however strong 
field (by the criterion r g /r). 
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The Newtonian limit (weak field conditions) is given by a linear approxima- 
tion of the exponential function 



{di/dOf = (1 - 2a r ) + 2a r £ - £ 2 



(51) 



where ay = t 9 /to7o/5o * s the a criterion in the relativistic case. 

Here, we recognize the Thomas precession proportional to r s /2r in the 
gravitational field due to the point mass M creating the potential energy 
field r g /r. In the strong field, the potential is 



Therefore, in the Newtonian limit, the relativistic Thomas precession is 
caused by the SR Kinematics. Next order corrections will be due to the 
deviation of the gravitational potential from the classical law V(r) ~ 1/r. 
Definitely, we do not confirm the GR precession proportional to 3r g /r . 
The theory predict a free radial fall without the deceleration 



Few more words about the theory. It requires a revision of the conven- 
tional relativistic concept of mass: the constancy of the proper mass in any 
field of forces. The revision results in an elimination of divergence of gravita- 
tional potential at r — > 0. It is shown that the introduction of the concept of 
variable proper mass also resolves the long standing problem of divergence in 
relativistic electrodynamics. For years, the problem has been tackled by the 
artificial procedure of mass/charge renormalization (by means of subtractions 
of infinite numbers). 

The theory adopts a field concept as an optical active medium for propa- 
gation of electromagnetic waves and material de Broglie waves. In this way, 
it has quantum connections needed for further development of relativistic 
gravitational field theory. The theory denies the GR "Black Hole" concept 
of gravitational collapse of matter and light. It explains all GR classical tests 
and has new practically verifiable predictions. 

8 Conclusion 

We raised questions, which are routine for any physical (small or great) find- 
ing: what and how it is observed, what and how it is predicted, finally, how 




(52) 




(53) 
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both sides are matched by a statistical method of empirical data treatment. 
The results of our investigation is summarized below. 

Astronomical observations of the GR perihelion advance. 

A direct observation of the effect is difficult because of its smallness on the 
background of classical perturbation and precession of the Mercury's orbit. 
The effect can be, in principle, assessed from the analysis of astronomical 
observations and the correspondingly adjusted ephemerides with the use of 
statistical methods. The ephemerides production is a multi-aspect complex 
work, and there had to be a spacial research program devoted to the Ein- 
stein's prediction of fundamental importance. In reality, there was no such a 
program. There are many published claims and statements, but there is no 
work which would present results of the whole problem study. 

Historically, sporadic initiative studies of individual researchers or groups 
were conducted. The goal mostly pursued was "to fit" the Einstein's effect 
value to the "anomaly gap" as close as possible. No rigorous methods of 
empirical data treatment from statistical theories were used by researchers 
for the effect verification. The implementation of the PPN formalism in 
the ephemerides systems made the situation worse because of uncontrolled 
unphysical shift in calculational results. 

Sadly, the opinion on the firm confirmation of the GR perihelion advance 
spread among the astronomical and physical communities. We state, how- 
ever, that the real precision of current ephemerides is unknown, and our 
knowledge about the effect verification in the astronomical observations re- 
mains, roughly, at the Le Verrier's level. 

GR prediction of the effect. 

The GR description of the perihelion advance is given in such a form that 
the effect cannot be discriminated from the huge fluctuating background of 
classically induced precession. Consequently, astronomers had to identify 
the effect as a difference of big fluctuating numbers in observed and cal- 
culated ephemerides values. As emphasized, they tried "to fit" theoretical 
ephemerides to the exactly predicted number. This methodology and the 
corresponding results cannot be termed an observational test. 

Rigor of the effect prediction is in a serious doubt. The Einstein's "ap- 
proximate" solution, when put back into the original equation, does not fit 
the equation to the precision better than the effect value. The GR-term, 
which is thought to be the cause of both the perihelion advance effect and 
the light bending effect, has no physical sense, as shown by Fock and in our 
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work. The Fock's work also shows that the equation of light propagation in 
a vicinity of massive object must have the linear (potential) term instead of 
the GR-term. The latter is of the next degree of smallness and should be 
neglected. 

Controversies in the GR theory arise in connection with gravitational 
properties of both the particle and the photon. The GR prediction of the 
particle deceleration in a radial free fall has no physical sense and most likely 
is wrong. The prediction of the bending of light is not valid for the same 
reason as in the perihelion advance case, while the predicted red-shift is 
inconsistent with the GR framework. These and other arguments raise the 
question about sufficiency and completeness of the GR physical foundations. 

Overall, we conclude that the claimed confirmation of the GR prediction 
of the relativistic perihelion advance is neither theoretically nor empirically 
substantiated. 
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A Rigorous Statistical Approach 

A.l Definitions and general formulation 

Information on the general statistical concepts and the corresponding method 
of measured data treatment are available in literature, for example [42] and 
elsewhere. A presented statistical method is applicable in general case of any 
active experiments, passive observations, and any kind of measurements, as- 
sessments of experiment informativeness, planing experiments, etc. Roughly, 
the problem is formulated as next. Suppose, some observable system is stud- 
ied. The empirical data include measured quantities and their statistical 
dispersions with possible correlations in the form of covariance matrices. 
One needs to assess system parameters in order to predict system character- 
istics, which are, in general, functions of spatial coordinates and time with 
constraints imposed on a set of parameters within the model. 

Definitions and denotations. 

X => Xi (i = 1, 2... J) - a set of model parameters (masses M, rrii and 
others); 

F(X) => F k (k = 1, 2... K) - a set of observables, including ephemerides; 
they are functions of time t. 

S{F\X) => S ki = [dF k /dXi] 

x=x ~ sensitivity coefficients. They are 
obtained in numerical calculations. 

D(X), D(F) => D rs - covariance matrices; they are quadratic symmetric 
positive-definite matrices that describe dispersions of parameters and their 
correlations before the experiment (prior information) and after (posterior 
information). They are obtained for the prior parameters D(Xq), the mea- 
sured quantities D(F) as well the corresponding posterior estimates D(X*) 
and D(F*). A size of D(X) is determined by a number of parameters /, and 
for D(F) the number of measured quantities K. 

Matrix denotations here are similar to that in the linear algebra, where a 
row and a column are considered vectors what gives a conventionally defined 
operations of their dot-product and matrix multiplication. The superscripts 
< T > and < — 1 > are used for transposed and inversed matrices. The 
method is designed to find the best estimates of model parameters' values 
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and their dispersions and correlations arising from empirical data taken into 
account. 

The starting point is the introduction of the so-called likelihood function 
V(F,X), the criterion of matching measured results F to the correspond- 
ing calculational quantities F(X). This is the normal distribution of the 
random quantity characterizing \0(X) — C(X)\ 2 in terms of sensitivity coef- 
ficients S(F\X) and the statistical weights (the inverse covariance matrices) 
for measured quantities D^ 1 (F) as well as the prior data D^iXo). The 
initial dispersions and correlations are found from covariance matrices D(F) 
and D(X ). 

Thus, one needs to solve the extremum problem for the statistical distri- 
bution function 

V(F,X) = Const exp [-(1/2)AF t D(F)- 1 AF] (54) 

where 

AF = [5F — S(F\X)5X] is a random deviation of quantities F(X) from 
observations, see next. 

5F = (F(X )—F) are differences between calculated and measured quan- 
tities (must be small deviations); 

5X = (X — X ) are random deviations of parameters X from their prior 
estimates (must be small deviations). 

Measured quantities are, in general, non-linear functions, or functionals 
(functions of functions) of parameters F(X). To be consistent with the law 
of normal distribution, a statistical procedure of parameter evaluation should 
be conducted with a system of linearized functions 

F k (X t ) = F k (X i0 ) + [dF k /dX t ] x=XM 5(X t - X l0 ) (55) 

or in matrix denotations 

F(X) = F(X ) + S(F\X) X=X JX (56) 

This is how the sensitivity coefficients come out. 

Prior data X are assumed to obey a normal distribution law too: 

V(X) = Const exp [-(1/2)5X T J D(X )" 1 5X] (57) 

where 8X = (X — Xq). The criterion of "smallness" means that deviations 
from "true" valued of the parameters due to uncertainties are small enough 
to validate a linear approximation of functionals F(X) at the initial point 
X = X . If the point was appreciably shifted during corrections, sensitivity 
coefficients had to be recalculated at a new point. 
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A. 2 Data statistical treatment in the Bayesian approach 



In the Bayesian approach to statistics, posterior and prior information is 
connected through conditional probabilities P(X|F)-7 :, (J r ) = V(F\X)-V(X). 
Hence, the posterior estimate is related to the distribution 



V{X\F) = Const V(F\X)-V(X) 



(58) 



or m our case 



V{X\F) = Const exp \-(l/2)5X T D(X )~ 1 5X 



(59) 



From the above formulas, the maximum likelihood principle leads to the 
posterior estimate of parameters 



X* = X + D(X ) SD(X )S 1 + D(F) (F - F(X ) 



-i 



(60) 



The introduction of prior information into the assessment procedure gives an 
opportunity to conduct data treatment continually by portions as soon as new 
observational data become available. The prior data (mean parameters and 
covariances) are those assessed prior to new results come into the data base. 
With new results taken into account, the prior data become the posterior 
ones till next new results come, and so forth. Comparing the statistical 
standard error with the difference (O — C) one can make conclusions on data 
consistency. 

To evaluate informativeness and statistical significance of measurements, 
one needs to examine the prior and posterior covariance matrices for sets of 
parameters as well as measured and predicted quantities. For parameters, it 

iS r i-l 

D(X*) = [^(Xo)- 1 + S T D(F)- 1 S] (61) 

It can be given in the equivalent form that can be used to avoid an inversion 
of large matrices 



D(X*) = D(X ) - D(X )S T SD(X )S T + D(F 



-i 



(62) 



The assessments of new quantities (functionals) Y at X = X with sensitiv- 
ities S(Y\X) are similarly conducted. 

To account for the rapid periodic and slow evolution of the system, one has 
to calculate sensitivity coefficients as functions of time. There is a problem 
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related to the coordinate system transforms, which lead to the change of 
sensitivity coefficients. They should be recalculated in a new coordinate 
system and properly used in statistical evaluation procedures. It looks that 
our sketch suggests a lot of work one could try to avoid without a loss of 
already achieved apparently great values. 

The last formulas are used in practical evaluations with the input data 
X , D(X ), F(X ), F(X), DF{X) and the output data X*, D(X*), F(X*), 
D(F*). Notice that the covariance matrices are used to assess the posterior 
parameters while the latter are not needed to assess the covariance matri- 
ces. The sensitivity coefficients are used at all stages of data treatment. In 
the discussed problem, an analytical model is not available in practice; the 
coefficients should be found in numerical calculations S(F\X) = AF/AX). 

It is more convenient to work with relative quantities, namely, the input 
x = (X/Xq — 1) and the output / = (F/F — 1), correspondingly. Then 
the sensitivity coefficients must be defined in the dimensionless form too 
S ki = i^/^ 1 x _ x i an d the output data (adjusted parameters) are x* = 
(X* /X — 1). All covariance matrices will be related to the corresponding 
dimensionless quantities. 

It should be noted that the sensitivity coefficients can be geometrically 
visualized in terms of the linear (Euclidean) ND vector space with the corre- 
sponding linear algebra used above. Each parameter P is associated with a 
unit vector in a basis of linear independent vectors so that a function (or 
a functional) Pfe(P) is characterized by a vector Ski being a linear superpo- 
sition of basis vectors. The covariance matrix describe a dispersion ellipsoid 
(standard errors and correlations) of elements in a subsets of vectors (mea- 
sured or predicted quantities, parameters, and others). A correlation of two 
quantities A and B are characterized by a cosine of angle between vectors 
A and B. Hence, A and B are 100 % correlated, if collinear. A decrease of 
the ellipsoid volume (tnhe determinant of a covariance matrix) in the output 
of the experimental data treatment characterizes the precision improvement. 
In this way, the concepts of informativeness and information are rigorously 
defined. We would like to put an emphasis on the Bayesian approach to the 
statistical treatment of empirical data, as opposed to the the "classical" least 
square method. The latter reduces the problem formulation to the solution 
of N algebraic equations yi = J2j s ij x j with N unknown quantities (param- 
eters) x n by inversion of the matrix of the system a => a^ 1 (the matrix of 
sensitivity coefficients). Often, one or several parameters happened to be "ill 
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determined" because measured data y are not sufficiently sensitive to them. 
This would result in the matrix determinant being close to zero, det(a) << 1 
(the unfolding problem). The method fails. The use of prior information 
on the parameters in the first place is a physically natural method of statis- 
tical regularization of the solution (it makes the problem formulation "well 
posed"). 

A. 3 Role of constraints and initial conditions 

Let us discuss the role of constraints imposed on a system. As a simple exam- 
ple, consider a nuclear reactor model described by a set of nuclear constants 
X having statistical dispersions D(X ) (in terms of a covariance matrix). 
Calculated characteristics F(Xq) of a reactor are subject to reassessment with 
improved precision as a result of experiments. In the list of main character- 
istics is the neutron multiplication (effective) factor K e f, calculated value 
of which depends on nuclear constants and parameters. In practice, the 
reactor is designed to make a neutron balance in a stationary state that 
is, K e f = Const. This requirement imposes a constraint on a variation of 
constants to be taken into account in the sensitivity coefficient calculations. 
Namely, a variable parameter q(X) (usually, a control rod position in the re- 
actor, a function of X) must be chosen to restore the condition the neutron 
balance when a single constant X; L is varied: 

= dhf/dXi + (dk ef /dq) ■ (dq/dX,) (63) 
The sensitivity coefficient with the restoration of K e f is 

S{F k \Xi) = dF k /dXi - (dk e f/dq) ■ (dq/dXi) (64) 

In the discussed ephemerides problem, the variable parameters of the plane- 
tary system are subjected to constraints supposed to be embedded into the 
model. They arise from conservation laws: the constancy of total energy E 
(the symmetry of time translation that is, the balance of potential and kinetic 
energy), the constancy of angular momentum L (the symmetry of 3-space 
isotropy that is, the absence of torque), and the constancy of linear momen- 
tum P (the inertial system property that is, the absence of acceleration of a 
test particle at the barycenter). They should be considered similarly to the 
above example of neutron balance conservation. 
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The existing practice of the ephemerides production suggests that every 
calculation of ephemerides should start with the "standard" multi-digit ini- 
tial values at a certain instant of time. From the statistical viewpoint, such 
a procedure makes a strong influence on the parameters' adjustment. The 
a priori data become subjected to the systematical off-set what results in 
diminishing of significance of new observational data. This fact partly ex- 
plains why results of calculations conducted in different ephemerides centers 
look stable and similar. It gives a false expression of a very high precision in 
descriptions of the Solar planetary system while the real precision remains 
unknown. 

The presented method suggests that the initial data should be used once 
at a starting time. After that, the deviations from the once set a priori 
database are calculated upon arrival of every new portion of observational 
data. A new precision is determined from the corresponding a posteriori 
covariance matrices. 

A. 4 Concluding remarks 

We described the method of observational data treatment based on a rigorous 
statistical theory. The reason for its presentation here is that in the existing 
practice of treatment of astronomical observation data the terms are used 
such as "precision" , "least square" , "fitting" , "sensitivity" and so forth, from 
the vacabluary of statistical theories. However, from what we lerned from 
literature related to the problem, those terms, as they are used in reality, 
have nothing to do with the statistical theory. 

So far, we consider the Newtonian iV-body model. An introduction of 
GR relativistic effects inevitably leads to complications of the problem for- 
mulation for, at least, two reasons: a) such quantities as potential and kinetic 
energies are out of the GR arsenal; b) the exact GR iV-body solution, which 
would have the Newtonian limit, does not exist. That is why we are critical 
of the PPN concept and oppose its introduction in the ephemerides systems. 
The basic relativistic effects are known without their GR connections and 
the corresponding corrections can be easily introduced. 

There must be a few levels of parameter categorization with regard to 
their influence. The first level: fundamental constants G, Co (in principle, 
they can be considered the exact numbers). The next level are basic physical 
characteristics such as point masses of the Sun and planets (the N body prob- 
lem). Further level would deal with most numerous and shaky parameters 
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responsible for the tidal forces due to body size, composition and rotation, 
space debris, as well. Finally, one can think about relativistic corrections of 
the model but no new parameters are needed (except for the PPN case). It 
is important for researchers to have an opportunity to handle the parameters 
by their levels (to be able "to turn" them "on" and "off", for example, to 
compare results of Newtonian model with the one including relativistic cor- 
rections). Because parameters are characterized by a variety of measuring 
units, it is very convenient to formulate the problem in a dimensionless form. 

The problem of initial conditions is emphasized in the context of ephemerides 
adjustment. In practical numerical calculations of ephemerides, several dif- 
ferent files of "standard" initial data with a high number of significant digits 
are used, which are special parts of the systems. When new observational 
data became available for further testing and correcting the ephemerides 
data base, the same initial data are often used. It is explained why, from the 
statistical point of view, such practice is inconsistent with the statistical ap- 
proach. Renewable prior data files must launch the system from any instant 
T. 

Our main reason for discussions of the role of statistical approach to 
astronomical observation treatment is the fact that the GR effect of the 
Mercury's perihelion advance is claimed to be successfully confirmed. In the 
present work, the results of our analysis of the problem show that the claimed 
confirmation is not true. Astronomers have to bear some responsibility for 
their practice of observation treatment. 
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